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Abstract 

We analyze the universal properties of a new two-dimensional quantum gravity 
model defined in terms of Locally Causal Dynamical Triangulations (LCDT). 
Measuring the Hausdorff and spectral dimensions of the dynamical geometrical 
ensemble, we hnd numerical evidence that the continuum limit of the model lies in 
a new universality class of two-dimensional quantum gravity theories, inequivalent 
to both Euclidean and Causal Dynamical Triangulations. 
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1 Causal quantum geometry in two dimensions 

The approach of Causal Dynamical Triangulations (CDT) [1] provides concrete 
evidence that one must include causal, Lorentzian properties in the nonpertur- 
bative gravitational path integral in order for the associated quantum gravity 
theory to possess a classical limit. This should be contrasted with purely Eu¬ 
clidean constructions where the path integral is taken over geometric conhgu- 
rations which represent four-dimensional “spacetimes”, but do not contain any 
information about time, light cones or causality. The problem with the latter 
appears to be that - quite independent of the elementary “building blocks” used 
for constructing individual Euclidean conhgurations - their nonperturbative sum 
is completely dominated by highly degenerate objects whose superposition never 
leads to a four-dimensional extended universe on macroscopic scales, no matter 
how one looks at it. 

The Euclidean precursor of CDT, so-called Dynamical Triangulations (DT), 
is a case in point. In both DT and CDT quantum gravity one looks for scaling 
limits of regularized path integral expressions, where curved geometries are rep¬ 
resented by triangulated, piecewise flat manifolds. However, the inhnite-volume 
limit of the Euclidean theory is dominated by degenerate phases with no obvious 
physical interpretation in terms of general relativity [2], and phase transitions are 
of hrst order 011]. By contrast, one of the phases of CDT quantum gravity is 
characterized by a quantum geometry which on large scales exhibits properties 
of a four-dimensional de Sitter space [5], and the phase space of CDT contains 
a whole line of second-order critical points [6], which are being investigated as 
natural candidates for defining the searched-for continuum theory [7|. 

At the inception of the CDT approach, the toy model version of the theory in 
two spacetime dimensions played an important role: the nonperturbative CDT 
path integral over geometries can in this case be evaluated analytically [8], and 
the resulting two-dimensional quantum gravity theory compared with the much- 
studied theory of two-dimensional Euclidean quantum gravity, which likewise can 
be formulated and solved exactly in terms of dynamical triangulations (see [9] for 
reviews). The two theories turn out to be inequivalent, and are characterized by 
different critical exponents (see [lO] for a comparison). In terms of continuum 
formulations, they he in the universality class of Liouville quantum gravity m for 
DT and two-dimensional, projectable Hofava-Lifshitz gravity [12] for CDT. These 
models provide the first completely explicit example that “signature matters” in 
the context of the nonperturbative gravitational path integral, aka the “sum over 
histories”. As mentioned above, we now have good evidence that the same is true 
for the physically relevant case of quantum gravity in four dimensions. 

This paper will expand on the theme of two-dimensional quantum gravity as 
an interesting testing ground for quantum gravity proper, where both analytical 


2 


and numerical solution methods can be employed and compared. We will study 
the two-dimensional implementation of a recently introduced version of CDT, 
which goes by the name of “Locally Causal Dynamical Triangulations (LCDT)” 
[T3] or “CDT without preferred foliation” H. The path integral in these CDT 
models is performed over a class of piecewise flat Lorentzian geometries that is 
enlarged compared to standard CDT quantum gravity. The geometries are still 
causal, in the sense of having a well-dehned light cone structure everywhere, but 
are not required to have the preferred (discrete) proper-time slicing character¬ 
istic of standard CDT conhgurations. An in-depth numerical investigation of 
locally causal DT in three spacetime dimensions found nontrivial evidence that 
key results of CDT quantum gravity, including the volume distribution of three- 
dimensional de Sitter space are reproduced in this generalized causal theory [H]. 

This is an important and concrete piece of evidence that for a judicious choice 
of the bare coupling constants of the theory, LCDT and CDT quantum gravity 
lead to equivalent continuum theories. We would like to investigate whether the 
same is true in two spacetime dimensions. Although this toy model is arguably 
even less representative of full gravity than the three-dimensional model, the 
properties of “quantum geometry” are much simpler to analyze in dimension two 
and may give us a hint of why the two causal theories are equivalent or not, as 
the case may be. 

Since in terms of its conhguration space the locally causal model lies in be¬ 
tween DT and CDT quantum gravity, solving it will give us a better understand¬ 
ing of the universality classes of theories of quantum geometry in two dimensions. 
The CDT universality class has so far proven to be quite robust: inclusion of a 
higher-curvature term na, a decoration by arches along spatial links (tantamount 
to including a restricted class of “baby universes” ) naiiz], an explicit inclusion of 
a hnite number of baby universes within a string held-theoretic setting based on 
CDT [IS] (see also [IS]), or starting from a conceptually rather different Hamil¬ 
tonian “string bit model” 1201 all lead to the same scaling limit. In the absence, 
to date, of an analytic solution of locally causal DT in two dimensions, we will 
present below the results of a numerical investigation. We have examined sev¬ 
eral observables and measured two critical exponents, the expectation values of 
the Hausdorff and the spectral dimension of quantum spacetime, to try to un¬ 
derstand whether they coincide with those of DT or CDT quantum gravity, or 
perhaps signal yet another, new universality class of two-dimensional quantum 
geometry. 

We begin our analysis by introducing the locally causal DT model and its 
geometry in Secs. and and outline the set-up for the Monte Carlo simulations 
in Sec. IH Closed timelike curves and their role in LCDT are described in Sec. [H 
Sec. deals with observables, including the volume prohle, a characterization of 
the quantum geometry in terms of minimal loops, and our results for the Hausdorff 
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and spectral dimensions. Onr conclnsions are presented in Sec. The Appendix 
contains some details on the geometric shape of typical LCDT confignrations. 

2 Locally causal dynamical triangulations 

We will hrst derive an expression for the gravitational action in locally cansal DT 
in two dimensions. Onr starting point is the two-dimensional Einstein-Hilbert 
action in the continnnm, 

S = Kj d^xy/\i\{R - 2A), (1) 

where k is the inverse of Newton’s constant, A the cosmological constant, and 
g denotes the determinant of the Lorentzian spacetime metric The integral 
over the scalar cnrvature R is topological and will not play a role in the path 
integral constrnction, since the spacetime topology will be held fixed. Dropping 
the i?-term leaves us with just the volume term. Absorbing the (dimensionless) 
gravitational constant into the cosmological constant A, the path integral and its 
regularized counterpart in terms of dynamical triangulations read 

Jvy] 

where the (formal) integration over diffeomorphism equivalence classes [g] of met¬ 
rics on the left-hand side has been replaced by a sum over inequivalent triangula¬ 
tions T with the usual DT measure involving the order C{T) of the automorphism 
group of T. The constant A on the right-hand side is the bare cosmological cou¬ 
pling of the regularized theory, and V 2 (T) denotes the spacetime volume of the 
triangulation T. 

As described in |13], LCDT in 1-1-1 dimensions uses two types of triangular 
Minkowskian building blocks, to allow for the construction of geometries without 
a preferred time foliation (see Fig. [^. The usual CDT simplex with one 
spacelike and two timelike edges is supplemented by another two-simplex Agst 
with one timelike and two spacelike edges. The squared edge length of all spacelike 
links is l‘l = a?‘ and that of timelike links 11 = —aa?, in terms of the lattice cutoff 
a and the ratio a > 0 of the two quantities. To determine the spacetime volume 
V 2 {T) of a triangulation T assembled from these building blocks, we simply need 
to count their numbers Nstt{T) and Nsst{T) and compute the volumes of both 
Astt and Agst- The latter are determined in a straightforward way from the values 
of the edge lengths of the two triangles, and are given by 

vo\{Agtt) = -^- a , vol(A^st) =- - - a . (3) 


^-iA J<fixy/\g\ 
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Figure 1: The two elementary building blocks of (l+l)-dimensional LCDT, with 
light cones indicated: stt-triangle Agtt with one space- and two timelike edges 
(left) and ssf-triangle Aggt with one time- and two spacelike edges (right). The 
colour-coding for spacelike edges is blue, and for timelike ones red. 



Figure 2: Two time orientations are possible for each triangle type, as indicated 
by the future-pointing arrow assignments. 

To be able to perform Monte Carlo simulations of the system, we analytically 
continue the parameter a to —a in the lower-half complex plane according to the 
usual CDT prescription [21]. The resulting real expression for the Wick-rotated 
regularized path integral in two dimensions is 



Note that for both triangle volumes to be positive, a must satisfy the inequality 
l/4<a<4. The limiting value a = 1/4 corresponds geometrically to a collapse of 
the stt-triangles to zero volume, whereas for a = 4 the sst-triangles are collapsed. 
In the isotropic case, q; = 1, both triangles after Wick rotation are equilateral and 
identical. All LCDT simulations presented below were performed for q; = 1. 

Building blocks of the two triangle types are assembled into simplicial mani¬ 
folds T by “gluing” them together pairwise along boundary edges, where a time¬ 
like edge can only be glued to another timelike edge, and a spacelike edge only to 
another spacelike edge. Note that with respect to some overall flow of time, each 
of the two building blocks can appear with two different time orientations, which 
can be indicated by arrow assignments as illustrated in Fig. By dehnition, all 
arrows are future-pointing. Note that once a single triangle in a triangulation has 
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Figure 3: Local causality implies the existence of one future and one past light 
cone at each point. The edges meeting at a vertex are always arranged in four 
groups of alternating type, time- or spacelike (with at least one edge in each 
group), as depicted. In the hgure, “time” is pointing upwards, and the thin lines 
indicate the light cones located at the central vertex. 

been given a specihc time orientation, the orientation of its neighbours, and of 
its neighbours’ neighbours, etc. is also hxed, because the arrows on shared edges 
have to match. 

Local causality is incorporated in the gluing rules by stipulating that before 
the analytic continuation there should be exactly one future and one past light 
cone at each interior point of the triangulation [H]. This condition is always 
satished at an interior point of a triangle, because up to diffeomorphisms the 
metric by construction is given by the Minkowski metric. It is also satished at 
points along edges where two triangles meet, unless the point happens to be a 
vertex, as can be seen by inspecting the geometry of the building blocks in Fig. 
1^ The requirement is only nontrivial at the vertices of the triangulation. When 
expressed in terms of the edges meeting at a vertex, it implies that they should 
come in four groups of alternating type (time- or spacelike) when going around 
the vertex once (Fig.|^, which imposes corresponding restrictions on the triangles 
meeting at the vertex. 

The generic vertex structure in ordinary CDT - which uses only building 
blocks of type lS.stt ^ is also of this type, but by construction there is only a single 
spacelike edge each on the left and the right of the light cone, and the pair of 
these spacelike links forms part of a preferred slice of constant integer time. It 
is precisely the generalized vertex structure in locally causal DT that allows for 
conhgurations without this preferred time slicing. Fig. [^illustrates the difference 
between a piece of causal DT and one of locally causal DT. 
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Figure 4: Two strips of dynamically triangulated spacetimes with an initial (bot¬ 
tom) and a final spatial boundary (top), (a) In standard CDT only stt-triangles 
are used. There are two spatial edges meeting at each internal vertex, giving rise 
to the characteristic preferred proper-time slicing in terms of consecutive lines of 
spatial links, (b) LCDT works with both stt- and sst-triangles, allowing for a 
more general vertex structure, without preferred slicing. 

3 Properties of LCDT in two dimensions 

Our task will be to investigate the properties of the path integral 0 in the 
continuum limit, where the sum is taken over an ensemble of triangulations of 
hxed topology, obeying local simplicial manifold condition^ and with the vertex 
structure of locally causal DT described above. At this stage, there is no known 
exact solution of the continuum dynamics of LCDT quantum gravity in two 
dimensions; one difficulty is precisely the absence of a distinguished notion of 
time in terms of the lattice structure itself, which prevents the straightforward 
introduction of a transfer matrix used previously in solving CDT [HI HSl [161 IE] • 
Since the configurations have two different kinds of edges, they can be thought 
of as a particular kind of two-coloured graphs, whose properties one may try to 
understand in a systematic way in the sense of enumerative combinatorics. The 
subgraph consisting of spacelike links only has the form of a stack of “bubbles”, 
which on their inside are decorated with timelike links (see (Haim for definitions 
and discussions of this substructure). The LCDT model may be supplemented 
by additional conditions, which restrict the type of (self-)overlaps among these 
bubbles that are allowed. Since the bubbles are extended structures, these addi¬ 
tional rules have a nonlocal character. They are motivated by the finding that 
for spatially compact boundary conditions in two dimensions the local condition 
of vertex causality described above does not imply global causality in the sense of 
the absence of (a specific class of) closed timelike curves (see Sec.j^for an explicit 
example and further discussion). In turn, the presence of such curves is related 
to the presence of overlapping bubbles. 

^in dimension 2: each internal edge is shared by exactly two triangles, and any two triangles 
share at most one edge 
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In our simulations both space and time will be compact. The topology of 
space will be a circle and time will be cyclically identihed, which means that 
spacetime is topologically a torus T^. Since we are primarily interested in “bulk” 
properties of the geometry, this choice is technically convenient: all vertices are 
interior vertices, on which vertex causality will be imposed, and the action is the 
one appearing in the path integral expression Q, without the need for adding 
any boundary terms. 

The functional form of the Euclidean action, schematically given by S' = 
ciNstt + C 2 Nsst, for two positive constants ci, C 2 , is the most general one linear in 
“counting variables”. These are the variables counting simplices of a particular 
dimension and type in a given triangulation T: the numbers Nq{T) of vertices, 
Ns(T) of spacelike edges and Nt(T) of timelike edges, as well as the numbers 
NsstiT) and NgttiT) already introduced earlier. Our statement follows from the 
fact that these hve variables are subject to three constraints, 

NQ — Ns — Ni,-\-Nsst'^Nstt=0, iVsst + 2iVsti —2iVi= 0, Nstt~\-‘2‘Nsst~‘^Nt= 0, (5) 

which must be satished on each conhguration T with torus topology. 

We hnally note that, at least in the absence of additional constraints on the 
bubble conhgurations, toroidal boundary conditions introduce a duality into the 
two-dimensional LCDT system. The duality transformation consists in swapping 
simultaneously the assignments “timelike” and “spacelike” of all edges in a given 
triangulation, which will convert all stt-triangles into sst-triangles and vice versa. 
An admissible triangulation (one that satishes the local conditions of a simplicial 
manifold and vertex causality) will under this transformation be mapped to an¬ 
other admissible triangulation with the same topology, with the roles of time and 
space interchanged. Of course, for a^l the Boltzmann weights of a triangulation 
and its dual will in general be different. 

4 Numerical set-up 

We have used Monte Carlo techniques to sample the partition function or Eu- 
clideanized path integral (|^ of locally causal dynamical triangulations, and com¬ 
pute expectation values of selected observables. An important ingredient are a 
set of Monte Carlo moves, which take the form of local changes in the simplicial 
geometry, and are designed to get us around the conhguration space of the model 
by way of a Markov process. The four types of move we have worked with will be 
described below. Further technical details on their implementation may be found 
in [23]. Additional references on Monte Carlo moves in the context of causal 
dynamical triangulations are 



Figure 5: Two types of Monte Carlo moves used in locally causal DT: (a) example 
of a (0,2)-move and its inverse, also called a “link collapse”; (b) example of a (2,2)- 
or flip move and its inverse. 


The first type of move is a generalization of the (0,2)-move used in CDT, in 
which two adjacent stt-triangles that share a spacelike link are created simul¬ 
taneously. The local starting conhguration consists of a pair of timelike links, 
belonging to opposite sectors of a vertex (one link from inside the past and the 
other from inside the future light cone), see Fig. [^a). This move is compatible 
with the time slicing of CDT geometries. Since this compatibility is no longer 
a requirement in locally causal DT, we will also use the colour-reversed counter¬ 
part of this move, where two adjacent ssf-triangles that share a timelike link are 
created from a pair of spacelike links, this time from opposite spacelike sectors of 
the light cone at the central vertex. 

Also the second type of move, the (2,2)- or “flip” move generalizes a local move 
already employed in CDT. It consists in flipping the diagonal inside a rhombus 
made of a pair of adjacent triangles. The version depicted in Fig. |^b) is the one 
also permitted in CDT. In our simulations of locally causal DT, we will in addition 
use the move with the opposite assignments of time- and spacelike edges. These 
are the only two flip moves compatible with vertex causality, if the character of 
the flipped edge remains unchanged. Two more flip moves are possible if the 
flipped diagonal link is allowed to change from time- to spacelike or vice versa. 

Another type of move we have used in the simulations is a (2,4)-move, where 
the starting point is again given by a pair of adjacent triangles. A new conhg¬ 
uration with identical boundary is obtained by “subdividing” the rhombus with 
another diagonal, thereby creating a four-valent vertex at the centre, see Fig.|^a) 
for an example. In order for vertex causality to be satished at the new vertex, 
the added diagonal has to be of opposite (time-/spacehke) type to the one al¬ 
ready present. Eight variations of the (2,4)-move (and its inverse) are possible, 
depending on the type and orientation of the initial triangle pair, but six of them 
are equivalent to performing a (0,2)-move, which was already discussed above. 

Lastly, we employ a “pinching” move, which is entirely new and not present 
in standard CDT, and was previously described in [13]. The initial local conhg- 
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Figure 6: The two remaining types of Monte Carlo moves used in locally causal 
DT: (a) example of a (2,4)-move and its inverse; (b) the pinching move and its 
inverse. 


uration for this move looks like a piece of regular CDT conhguration, consisting 
of four stt-triangles forming a strip bounded above and below by line segments 
made out of spacelike links. The move pinches those two segments together in a 
single point, resulting in a pair of sst-triangles, as illustrated by Fig. [^b). In the 
simulations we also use the colour-reversed version of this move. 

In all cases, it is understood that whenever one of these moves is proposed 
by the computer algorithm, it will always be rejected if it violates either vertex 
causality or the simplicial manifold condition. The (overcomplete) set of these 
moves is likely to be ergodic, but we do not have a formal proof at this stage. 
The explicit proof may depend in subtle ways on the details of the ensemble, in 
particular, on excluding classes of bubble conhgurations associated with specihc 
closed timelike curves that lead to unwanted global acausal behaviour. 

We have run several kinds of cross check on the Monte Carlo simulations: 
hrstly, that the acceptance rates of moves and their inverses are approximately 
the same, and secondly, that the frequency of occurrence for conhgurations with 
very small volume is compatible with the frequency predicted by the Boltzmann 
distribution. Lastly, our set-up has an easy way to implement a CDT limit, which 
we can use to cross-check measured CDT results for the dynamical dimensions 
with the theoretical results available. 

Note that apart from the hip move, all Monte Carlo moves described above 
alter the number of triangles in the triangulation, and therefore its two-volume. 
Since it is convenient from the numerical point of view to keep the total volume 
hxed, at least approximately, we use the standard DT prescription where the 
volume V 2 is allowed to vary in a narrow interval around a hxed target volume 
¥ 2 ^^. This is achieved by adding a quadratic term 6 (14^°—V^)^ to the action, with a 
parameter 5 > 0 determining the width of the interval]^ We tune the cosmological 
constant A such that the measured volumes are distributed symmetrically around 

^Note that for our standard choice a=l, V 2 is proportional to the number N 2 of triangles. 
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Figure 7: A piece of LCDT, with an initial and a final spatial boundary (blue 
edges, bottom and top). When the two spatial boundaries are identified as indi¬ 
cated by the letters, the future-oriented timelike curve running from bottom to 
top (timelike edges with arrows) becomes a timelike cycle, as defined in the text. 

V 2 ^\ and only collect data from triangulations that have precisely this target 
volume. The simulations are run such that there is about one sweep of length 
L K. 10® between successive measurements. 

As mentioned before, there is a straightforward way to obtain ordinary CDT 
simulations in our set-up. It consists in setting a = 1/4 in the action 0 , while 
maintaining the constraint N 2 = Ngst + Ngu = constant. Since the term pro¬ 
portional to Ngtt in the action now vanishes, stt-triangles can be created “at no 
cost” during the Monte Carlo simulation, while the number of ssf-triangles will 
diminish accordingly, thereby lowering the value of the Euclidean action. As a 
result, the triangulations quickly become pure CDT configurations, consisting 
only of building blocks Agu- Below in Sec. [^we measure the Hausdorff and spec¬ 
tral dimensions of CDT quantum gravity. In addition to setting a = 1/4, we 
will disable the ssf-triangles completely, to make sure that any fluctuations with 
nonvanishing Nggt are eliminated. 


5 Closed timelike curves 

For most of our measurements, the ensemble of locally causal triangulated ge¬ 
ometries on which the dynamics takes place will consist of spacetimes of torus 
topology satisfying simplicial manifold conditions and vertex causality. However, 
for some purposes, when considering the time evolution of observables, it is con¬ 
venient to restrict this ensemble further, because of the appearance of a particular 
class of closed timelike curves. 

To explain this issue further, we introduce the notion of a timelike cycle in a 
locally causal triangulation. By this we shall mean a contiguous set of timelike 
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Figure 8: A piece of locally causal DT, with an initial and a final spatial boundary 
(blue edges, bottom and top). Compactifying space as indicated by the letters 
results in a closed, future-oriented timelike curve (timelike edges with arrows), 
which intersects neither the initial nor the hnal boundary. 

links which together form a non-contractible loop of topology S^, without self¬ 
crossings or self-overlaps (Fig. [^. In addition, whenever the loop passes through 
a vertex, the two timelike links of the loop meeting at the vertex must lie in 
opposite light cones, never inside the same (half of the) light cone. A spacelike 
cycle is dehned analogously in terms of spacelike edges, and it is also required to 
be non-contractible. In the spacelike case, the cycle has to cross at each vertex 
from one spacelike sector outside the light cone to the opposite one. 

In standard CDT in H-l dimensions with its preferred time slicing (c.f. Fig. 
i>). spacelike cycles only exist when space is compactihed to a circle. In this case 
they simply coincide with the one-dimensional spatial slices at integer proper 
time. Likewise, timelike cycles only exist in CDT when time is compactihed, in 
which case they can be thought of as a particulai0 lattice realization of closed 
timelike curves. As a consequence of how they traverse the light cones at vertices, 
they are also time-oriented, either in positive or negative time direction. As 
usual, the reason for compactifying time in CDT simulations is merely one of 
convenience, and the appearance of closed timelike curves is an inevitable side 
effect, which is not expected to have much inhuence on the measurement of “bulk” 
observables like the dynamical dimensions considered belowj^ 

The closed timelike curves we are primarily concerned about in LCDT are not 
the ones winding around the compactihed time direction, but around the spatial 

^They are particular in the sense that one could also consider paths that run not only along 
the edges of a triangulation, but also through the interiors of triangles. 

"^Obviously, in any concrete implementation the dependence of observables on boundary con¬ 
ditions and other finite-size effects should always be monitored. Although CTCs are unphysical 
classically, their inclusion in the regularized path integral does not a priori imply unphysical 
behaviour of the final continuum theory. 


12 








direction, and which would still be present if spacetime was a cylinder I x 
(with compact spatial slices) instead of a torus, see Fig. for an example. In 
the remainder of this work, when we talk about closed timelike curves (CTCs), 
we will mean only this class of timelike cycles. It turns out that local vertex 
causality does not preclude the presence of these curves, although it does prevent 
the occurrence of contractible timelike loops [22] • One way of hnding them is by 
running an algorithm that assigns time labels to vertices of a given locally causal 
DT. Of course, an explicit choice of time has to be made in the LCDT model, 
because - unlike in usual CDT - there is no preferred lattice substructure one 
can refer to as a natural time label. 

The prescription for assigning time labels to vertices in a given geometry we 
have used in the present work is to pick a spacelike cycle in the geometry and 
dehne it to be “space at time t = 0”. Vertices lying in the future of this slice 
are then successively assigned time labels, which are computed as the average 
distance of the vertex v to the initial slice along any oriented timelike path from 
the slice to v (see [23] for more details on the algorithm). The distance along 
any given path is simply given by the number of timelike links it contains, the 
so-called link distance. Note that the time label of a vertex will in general not be 
an integer. Once the vertices are labelled, one can in a straightforward way also 
associate time labels with edges or more extended regions like spacelike cycles 
by averaging over the time labels of the vertices contained in them. In standard 
CDT, this prescription reproduces the usual integer proper-time slicing. 

The algorithm implementing the vertex labelling breaks down when it en¬ 
counters a CTC, like that depicted in Fig. because for any vertex lying on the 
curve or in its future there will be inhnitely many timelike paths connecting it to 
the initial slice. We conclude that in LCDT, at least when space is compactihed, 
local causality does not imply global causality in the sense of an absence of CTCs. 

Whether or not the presence of closed timelike curves has any consequences 
for the continuum limit of the model is a priori unclear. We have performed a 
number of measurements to get a better quantitative idea of how many CTCs 
there are, depending on the size of the triangulation. Fig. shows two his¬ 
tograms of the frequencies of disjoint CTCs (CTCs without mutual overlap]^. 
We observe that the typical number of disjoint closed timelike curves present in a 
given conhguration goes down signihcantly when the volume is increased tenfold 
from W = 10.000 to N 2 = 100.000 triangles. On the other hand, the ratio of 
triangulations which contain any CTCs at all increases from about 21% to 30%. 
We therefore have no indication that CTCs disappear completely as the volume 
grows. One contributing factor is presumably that it takes much longer to break 
up a CTC by a local Monte Carlo move when the CTCs become highly diluted 

® These curves form a subclass of CTCs; the number of all CTCs can be considerably larger, 
especially when there are many disjoint CTCs. 
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Figure 9: The frequency of disjoint timelike loops in 48.000 sweeps. The upper 
histogram is for A ^2 = 10.000 and the lower one for iV 2 = 100.000. 


in a triangulation. 

6 Observables 

Having introduced both the theoretical and numerical set-up of locally causal 
dynamical triangulations in 1-1-1 dimensions, we will now discuss the measure¬ 
ments of several observables in this model of quantum geometry, concentrating 
on the isotropic case a = l. Our main aim is to understand whether the model’s 
continuum limit coincides with that of either Euclidean or causal dynamical tri¬ 
angulations in two dimensions. It would be exciting if LCDT constituted a new 
universality class, but this seems a priori less likely because of the apparent 
scarcity of universality classes among two-dimensional models of pure geometry 
without matter coupling. 

If LCDT quantum gravity were to lie in one of the two known universality 
of DT models in two dimensions, our best guess at this stage would be that it 
is equivalent to CDT, for two reasons: hrst, there is good evidence that this is 
true in three spacetime dimensions [H], in the sense that there one hnds in both 
models a phase whose ground state has the scaling properties of a Euclidean de 
Sitter universe. (Of course, this by no means constitutes a proof that the same 
happens in two dimensions, which differs in terms of both its geometric degrees of 
freedom and its phase structure.) Second, the difference between DT and CDT in 
two dimensions has so far been explained in terms of the absence of baby universes 
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in the latter uni [21113 Since the condition of vertex causality in locally causal 
DT suppresses the light cone degeneracies characteristic of topology change and 
therefore of the creation of baby universes, the LCDT model seems closer to CDT 
than to DT, where baby universes dominate. 

6.1 Volume profile 

We begin by examining the so-called volume prohle of a typical geometric conhg- 
uration generated by the Monte Carlo simulation of LCDT. The volume prohle 
is simply given by the size of the spatial volume as a function of time. Since 
by construction the model does not have a distinguished notion of time in terms 
of some lattice substructure, we will make use of the notion of time introduced 
in Sec. above. For this purpose, a spatial slice at hxed time is any spacelike 
cycle - as dehned at the beginning of Sec. and its time label is obtained by 
averaging over the time labels of all of its vertices. The volume of a spatial slice 
is the number of links contained in it. To determine the complete volume prohle 
of a spacetime conhguration, one has to identify all of its spatial cycles, which 
is a nontrivial task. Recall that unlike in CDT, in LCDT diherent spatial cycles 
can cross and overlap along some subset of spatial edges. 
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Figure 10: Sampled LCDT volume prohle for W = 10.000, where t denotes the 
time label of a spacelike cycle and the discrete spatial extension is the number of 
spacelike links in the cycle at a given t. 

®More precisely, it is the absence of the possibility for baby universes to proliferate without 
limit; a limited, controlled presence of baby universes is compatible with two-dimensional CDT, 
as demonstrated by the model of generalized CDT [181 HI] • 
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For simplicity, we have considered sampled volume profiles instead of complete 
ones; for a geometry of volume N 2 = 10.000 we have randomly sampled 100 spatial 
slices, and for each slice determined its volume and time label. Fig. [To] shows one 
such sample, to illustrate the situation. We note that some time labels are very 
close to each other, which indicates that they probably share one or more spacelike 
links. The sample shows large volume fluctuations within small time intervals, 
without any discernible overall shape. This hnding is only qualitative, but it is 
compatible with the typical, strongly oscillating volume profiles encountered in 
simulations of two-dimensional CDT [25]. 

6.2 Behaviour of minimal loops 

A phenomenon that will potentially affect the measurement of dimensions dis¬ 
cussed below is the overall shape of the toroidal conhgurations. We saw in the 
previous subsection that their volume prohles seem to be strongly fluctuating, 
and are comparable to what one hnds in two-dimensional CDT quantum gravity. 
Large fluctuations are commonplace in two dimensions, because there is only a 
single coupling constant (the cosmological constant) which sets the scale of both 
the spatial “universe” and its quantum fluctuations 00 This is in line with the 
fact that general relativity in two dimensions is trivial. 

A new feature in LCDT is the variable length of the conhgurations in the time 
direction, since by construction the hxed time slicing is absent. As a result, both 
the spatial extension of the universe and its time extension - determined by the 
prescription used for measuring the volume prohle, say - will evolve dynamically. 
For example, the torus may become very thin in one of its directions, an ef¬ 
fect which may be quantihed by monitoring the length of closed non-contractible 
curves. While in CDT individual slices of constant time can become very short 
(the minimal length of a spatial compatible with manifold conditions is at¬ 
tained by cycles of three links), the probability for this to happen can be made 
very small by choosing the total time extent Itot and the total volume N 2 suit¬ 
ably. By contrast, in LCDT it can in principle happen that the torus becomes 
uniformly thin in one of its directions, even for large N 2 , if this is dynamically 
preferred. 

The relevance of this for the measurement of dynamical dimensions is the 
possible appearance of hnite-size effects, even when the total volume is large. For 
example, this happens when the paths of random walkers - used to determine the 
spectral dimension - start winding around the torus more than once. To obtain 
an estimate of the size of this effect we have set up an algorithm which searches 
the tori for minimal non-contractible loops. It does not distinguish between time- 

^The unique length scale of the quantum theory is where A is the dimensionful, 

renormalized cosmological constant. 
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Figure 11: Evolution of the shortest and longest loop lengths from a sample 
of minimal non-contractible loops through randomly chosen vertices on a given 
triangulation, as a function of Monte Carlo time t, and at discrete volume N 2 = 
100.000. Largest t-values are just before the onset of thermalization. 


and spacelike links, which implies that the minimal loops found can be made up 
of any link types. 

The algorithm consists of the following steps: on a given LCDT conhguration 
T, pick a vertex v. Starting at v, perform a breadth-first search until the area 
searched starts to self-overlap at some other vertex v'. Next, determine whether 
the minimal closed curve c through v and v' obtained in this way is contractible 
or not, by starting breadth-first searches on either side of c. If those searches start 
overlapping, c is a non-contractible closed curve through v of minimal length and 
we record its length. Because this procedure is rather costly in computational 
terms, we do not repeat it for every vertex of T, but only for a sample of 1.000 
randomly chosen initial vertices on T. For a given triangulation T we there¬ 
fore end up with a sample of locally minimal loops (“locally” because they pass 
through prescribed vertices on the torus). 


To obtain Fig. IT, we have performed the sampling of minimal loops at each 
step during the thermalization of a Monte Carlo simulation of 280 time steps, 
corresponding to 280 sweeps, for a triangulation size iV 2 = 100.000. At each step, 
we plot only the shortest and longest minimal loop length of the sample; all other 
minimal loop lengths lie in between these two values. Since our samples are quite 
large, the lower curve is probably a good indicator of the global minimum of the 
length of non-contractible loops on T, which roughly speaking lies in the range 
20-35. This is far away from the kinematically allowed minimum of 3 and shows 
that the torus does not become very thin in some places]^ In addition, the fact 


®Note that our algorithm does not determine which of the torus directions any particular 


17 





that there is a signihcant distance between the upper and lower curves shows 
that the torus does not degenerate by becoming very long in one direction and 
uniformly short in the other. On the other hand, a triangle count of N 2 = 100.000 
according to (4) at a = l corresponds to a volume V = a/ 3A^2/4~43.300, where we 
have set the lattice constant to unity, a = l. Assuming the two torus directions 
are approximately of equal length, this corresponds to an average linear extension 
of the torus in the range of 150-200, of which the shortest loop length therefore 
is only a small fraction. It implies that fluctuations in the linear extension of 
the LCDT conhgurations are large, even if the total volume is also large. We 
will comment further on this characteristic feature of two-dimensional quantum 
gravity in Sec. 6.4 below. 

The overall conclusion is that at system size N 2 = 100.000 hnite-size effects for 
observables involving shortest (geodesic) distances should not play a role at least 
up to link distances of about 30, and for observables involving closed random 
walks at least up to about 500 stepsj^ Looking ahead to the measurement of 
dimensions presented below, this still leaves plenty of room for hnite-size effects 
on larger scales, but they are not quantihed easily above the thresholds just 
mentioned. 


6.3 Spectral dimension 

Dynamical dimensions, like the spectral and Hausdorff dimension, are important 
and popular examples of observables in models of nonperturbative quantum grav¬ 
ity because of their computational accessibility in many different contexts. A key 
insight is that the values of these dimensions do not have to coincide with the 
dimensionality of the triangular building blocks used to construct the regular¬ 
ized model if one takes a nontrivial, inhnite continuum limit, as we are doing. 
Furthermore, a familiar feature from studying graphs and fractals, namely, the 
fact that there exist “spaces” of non-integer dimension, is also encountered in 
systems of dynamical triangulations. This is not necessarily inconsistent from 
a physical point of view as long as the anomalous values of the dimensions are 
conhned to a highly quantum-fluctuating, non-semiclassical regime, typically at 
the Planck scale. Since there is no nontrivial classical theory of two-dimensional 
general relativity whose solutions might be recovered from a corresponding quan¬ 
tum theory in the limit as h —)■ 0, there are no a priori physicality constraints on 

minimal loop winds around. Also in this respect the information is distinct from that contained 
in volume profiles of the kind shown in Fig. [T0| 

®Note that we are not making a distinction between link distance on the triangulated lattice 
and link distance on the dual lattice. We have not determined the relative scale between these 
two notions of geodesic distance on typical LCDT configurations. The numbers given in the 
text should therefore be treated only as rough estimates. 
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the Hausdorff and spectral dimension of an ensembl^^ of DT configurations in 
two dimensions, causal or otherwise. 

For Euclidean DT in two dimensions, from theoretical scaling arguments the 
spectral dimension is 2 and the Hausdorff dimension 4 [26], which has also been 
corroborated numerically [27]. Invoking an equivalence between CDT conhgura- 
tions and tree graphs, CDT in 1+1 dimensions can be shown to have a spectral 
dimension of at most 2 and a Hausdorff dimension of almost surely 2 [2H] , the lat¬ 
ter in agreement with earlier theoretical mm and numerical [25] results. In the 
context of LCDT, we will hrst investigate the spectral dimension. In the section 
following this one, we will examine the Hausdorff dimension, which appears to 
be the quantity best suited to discriminating between the different universality 
classes. 

The hrst step in measuring the spectral dimension is to dehne a discrete 
diffusion process on a two-dimensional Euclideanized locally causal triangulation 
T. This takes the form of a random walk moving in steps of unit distance between 
the centres of neighbouring triangles as function of a discrete external diffusion 
time a, analogous to what was done in CDT in four dimensions pUlIHUj . In other 
words, the diffusion takes place along the edges of the trivalent lattice dual to T. 
Calling KTii, +; cr) the probability to go from triangle io to triangle i in a steps, 
satisfying +; f) = 1, fhe discrete diffusion equation on the triangulation 

T reads 

KT{i,io;a+ 1) = {1-x)KTii,io;a) + ^ ^ KT{j,io;a), ( 6 ) 

j n.n. of i 


subject to the initial condition Krii, io; (J = 0) = 6i^ig. The sum on the right-hand 
side of ([^ is over the three nearest neighbours j of triangle i, and y G [0,1] is a 
diffusion constant which allows for a non-vanishing probability (1 — y) that the 
random walker remains at the same triangle during a diffusion step. It is included 
merely for convenience, to somewhat smoothen out the discretization artefacts 
for short diffusion paths. In particular, there is an asymmetry between paths of 
even and odd numbers of steps (c.f. the discussion in [29]), with a corresponding 
oscillatory behaviour in the curve for dg that is also present in our Figs. [T^ and 
13 when one zooms into the region below a ~ 50. A diffusion constant y < 1 has 
been used previously when studying the spectral dimension in three-dimensional 
CDT [3l]. In our simulations, we have worked with y = 0.8 throughout. 

To extract the spectral dimension, we consider closed random walks, beginning 
and ending at a specihed triangle +• They enter into the calculation of the average 


^°When talking about “the Hausdorff dimension”, say, of (C)DT, we mean of course the 
expectation value of this quantity, measured in the ground state of the relevant ensemble. 
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return probability 


(7) 




1 

N2{T) 


KT{io,io;cr) 

io6T 


for a given triangulation T. The spectral dimension dg is obtained from the 
expectation value {P{a))N 2 of the observable Q in the ensemble of triangulations 
of hxed volume N 2 according to 


dg{a) 


o dln(P(a)) 

diner 


( 8 ) 


In practice, we perform a “double” random sampling, where for each randomly 
chosen triangulation we pick 10 times a triangle randomly as starting point zq 
for a random walk, and then repeat the process for at least 400 triangulations. 
For a diffusion process on classical, flat IR'^, the formula analogous to simply 
reproduces the topological dimension d, independent of a, but in the quantum 
context the behaviour of dg can be more complicated and, generally speaking, 
cr-dependent. By not denoting the A^" 2 -dependence in ([^ explicitly we mean to 
indicate that dg is determined in the limit of large volumes where this dependence 
gradually disappears. 

When the total spacetime is compact, the spectral dimension will always go 
to zero for sufficiently large a. This is a hnite-size effect which occurs when 
generic random walks become sufficiently long to wrap around space one or more 
times. We are primarily interested in the a-regime below this range. If the system 
develops a stable plateau below the scale where signiheant hnite-size effects kick 
in, we will refer to this constant value of dg as the spectral dimension of the 
underlying “quantum spacetime”. Note that if the system size is too small, a 
plateau will never form due to a dominance of hnite-size ehects. 

Since there are no published numerical results on the spectral dimension of 
two-dimensional CDT, and since it will be useful to have a point of reference 
for the measurements in LCDT, we will hrst present our results for the spectral 
dimension of regular CDT quantum gravity on a two-torus. As explained earlier, 
the reduction to pure CDT conhgurations is achieved by setting a = 1/4 in the 
action. To understand better the ehects of the discretization, we have used two 
diherent discrete versions of the dehning formula ([^ for the spectral dimension. 
Having determined the expectation value {P{a)) for integer a from the data, we 
have employed two diherent implementations in terms of hnite diherences. The 
standard choice is 


(1) ^ _ ln(P(a + l)) -ln(P(a)) 
ln(a + l)-ln(T 


(9) 


In addition, we have used the alternative form 



-2a 


f{Pia^_\ 

V (Pi^)) ) ’ 


( 10 ) 
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Figure 12: The spectral dimension of CDT as a function of the diffusion time 
a, measured at volume iV 2 = 100.000, and time extension f'roT = 80. The upper, 
green line is the dimension of eq. (9) and the lower, blue line is the dimension 
ds of eq. (10). Statistical error bars are too small to be displayed. 


which has the same continuum limit and the advantage that no expensive func¬ 
tions are required. 

The CDT results for the spectral dimension dg are displayed in Fig. 1^, for 
data taken at volume N 2 = 100.000 and time extension troT = 80. For a < 700 
there is a small discrepancy between the curves corresponding to the two different 
discretizations, giving us an estimate of the systematic error of determining dg for 
small values of a. For larger a, both curves merge into what is essentially a single 
plateau. The spectral dimension of CDT extracted from data on the plateau is 
dg = 2.02 ± 0.02, in good agreement with the expected value of 2. 

The curves for the spectral dimension for LCDT are shown in Fig. 


13 Qual¬ 


itatively the plot is similar to that of CDT, but the plateau is reached only for 
somewhat larger diffusion times a > 1.000. The numerical result for the spec¬ 
tral dimension is dg = 1.99 ± 0.02, which we regard as a convincing confirmation 
that the spectral dimension of locally causal DT is 2, like that for DT and CDT 
quantum gravity. 


6.4 Hausdorff dimension 

The Hausdorff dimension dh is a key quantity to discriminate between distinct 
universality classes of two-dimensional DT quantum gravity. The general idea 
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Figure 13: The spectral dimension of locally causal DT (with q; = 1) as a function 
of the diffusion time a, measured at volume N 2 = 100.000. The upper, green line 
is the dimension of eq. (9) and the lower, blue line is the dimension of 
eq. (10). Statistical error bars are too small to be displayed. 


is to relate the volume V of compact, connected regions in space - typically 
discrete analogues of geodesic balls around a chosen point - to their linear size, 
e.g. the radius r of the region, and to extract the leading scaling behaviour 
from (V^(r)) ~ For our purposes, we will use a “differential” version of 

this relation, where one monitors the one-dimensional volumes of spherical shells 
around a given triangle io or, equivalently, the number of (dual) vertices at radial 
distance r from a vertex of the lattice dual to a given LCDT configuration T. 
We define nir^io) as the number of triangles found at geodesic distance r from 
* 0 , where geodesic distance is defined as the (integer) length of the shortest path 
along edges of the dual lattice. We have ?7,(0,zo) = 1 and n(l,fo) = 3 for all zq, 
because each triangle has exactly three neighbours and therefore the dual lattice 
is trivalent. Every triangle of T will appear in exactly one of the shells, implying 
that = N 2 {T). The identihcation of the shells can be implemented 

as a modified breadth-first search, which keeps track of when a change of shells 
occurs. Averaging over all initial triangles io G T, we obtain the average shell 
volumes at radius r. 


n{r) 


N2 




*0 = 1 


( 11 ) 


In what follows, we will refer to the function n{r) as the shape of a triangulation. 
To extract the Hausdorff dimension, we have applied finite-size scaling to the 
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expectation value {n{r)) of the shape function (11 )j^ The simulations consisted 
of 48.000 sweeps each at volumes N 2 = 100.000, 200.000, 300.000 and 400.000, 
and were done for LCDT (at a = 1) and, for reference and comparison, also 
for CDT (corresponding to a = 1/4). The scaling ansatz for the radius and shell 
volume isr ^ X = N 2 and (n(r)) —)■ N 2 respectively. Individual 

data points (n(r)) were transformed into curves via spline interpolation and a 
Levenberg-Marquardt least-square fit was used to align the shapes 



Figure 14: CDT quantum gravity: Fit for best overlap of the rescaled shapes 
{n{r))/N 2 as function of the rescaled distance x = r/N 2 ^'^'", for Hausdorff 

dimension dh = 2.2. The extension in time direction was set to tTOT = 80. 


In maximizing the overlap we have taken into account all x-values where 
(n(r)) has at least half of its maximal value. We have also measured the “short- 
distance” Hausdorff dimension for small x by optimizing the overlap of the initial 
rising slopes of the curves (n(r)). In principle this dimension need not coincide 
with the global Hausdorff dimension dh we have been considering [32], but in 
our case there was little difference. Our results for the best overlap of the shape 
functions for the CDT case are shown in Fig. they correspond to a Hausdorff 
dimension dh = 2.2 ± 0.2, which is compatible with the known value of 2|^ Our 


11 


l/dn 


We also tried to extract the Hausdorff dimension from the scaling relation r{N 2 ) At 
for the average linear extension f defined in eq. (13), but this did not yield meaningful results 
because of the convergence issues to be described in more detail below. 

^^For CDT in two dimensions, the Hausdorff dimension has been measured previously [25j . 
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error bars are rather large, because the best £t depends quite sensitively on the 
x-range for which the overlap is optimized. 



Figure 15: LCDT quantum gravity: Fit for best overlap of the rescaled shapes 
{n{r ))/as function of the rescaled distance x = for Hausdorff 

dimension dh = 2.7. 


The analogous data for LCDT are displayed in Fig. 15 Maximal overlap 


is achieved for a Hausdorff dimension dh = 2.71 ± 0.2, which is far away from 
our conjectured CDT value of 2, and even further away from the DT value of 
4. We conclude that LCDT very likely lies not in the same universality class as 
DT. Contrary to our expectation, equivalence of LCDT and CDT appears to be 
excluded too. Instead, our measurements point towards LCDT lying in a new 
universality class, not hitherto seen in quantum models of two-dimensional pure 
gravity. This would be a truly interesting result, and it warrants another critical 
look at the strength of our evidence. 

As is apparent from Figs. M and lA, the quality of the overlaps is not very 
good. Could there be systematic sources of error that affect our results to the 
extent that they ultimately are not in contradiction with dh = 2 for LCDT? In 
other words, may we be underestimating our error bars signihcantly? It may be 
worth recalling that it took some time to nail down the Hausdorff dimension of 
two-dimensional DT quantum gravity numerically. In the words of the authors of 


with good results, from finite-size scaling of the distribution of spatial volumes. Because of the 
absence of a pre-defined time function, this method is not computationally feasible for LCDT. 
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|33j . early simulation results were “remarkably inconclusive” (see |23] for further 
references). The same work also used hnite-size scaling of the shape function 
to determine dh, with a £t quality somewhat similar to ours. Especially when 
using the dual lattice - as we are also doing in the present work - the Hausdorff 
dimension extracted this way was significantly off the mark [dh = 3.15 instead 
of the known, correct value 4). Of course, one should keep in mind that these 
simulations were performed for a geometric ensemble different from LCDT and 
for moderate lattice sizes N 2 < 32.000 only. On the other hand, the causal 
gluing rules of LCDT introduce a local “stiffness” in the conhgurations compared 
to DT, which is likely to require larger volumes to achieve numerical results of 
comparable quality. 

In the case of DT simulations, signihcant progress with respect to the conver¬ 
gence of hts was obtained by introducing a “phenomenologically fudged” scaling 
relation for the geodesic distance 123, namely. 


r -|- a 

+ b 


( 12 ) 


where a and b are two parameters meant to compensate lattice artifacts at short 


distances. We have also tested relation (12), but found that nonvanishing values 


for a and b steer the CDT results even further away from 2 and also increase the 
sensitivity to the choice of htting region. 

A generic feature of two-dimensional quantum gravity illustrated by the nu¬ 
merical difhculties already mentioned is the fact that in two dimensions quantum 
fluctuations are always large, even for large lattice volumes. This is different from 
CDT in higher dimensions, say, where the dynamics is governed by two scales: 
one macroscopic, related to the overall size of the universe, and another one mi¬ 
croscopic, setting the scale of quantum fluctuations. In two dimensions, there 
are no nontrivial classical solutions, and there is only a single scale, that of the 
quantum fluctuations. 

In this situation, it is therefore natural for hnite-size effects to generically 
be large, especially when there are non-contractible directions along which space 
can become “small”, as can happen for the torus topology used for LCDTp^ 
This is certainly relevant when measuring the Hausdorff dimension; when the 
geodesic balls centred at some triangle zq start wrapping around one of the torus 
directions, the interpretation of the scaling relation by which we extract dh will be 
affected, in the sense that only triangles not visited previously will be counted as 
belonging to a given radial shell. Of course, one can keep extracting the Hausdorff 
dimension regardless, but should be aware that it contains also global, topological 
information. 


^^Note that the DT simulations mentioned above use the topology of a two-sphere. Also this 
indicates the need to go to larger volumes in the LCDT case. 
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Figure 16: Development in Monte Carlo time t of the average linear extension 
f of a LCDT configuration T of volume N 2 = 50.000. Isolated peaks in f keep 
occurring, even as the number of sweeps becomes very large. (Note that the 
y-axis has an offset of 70.) 


To get further insights into the origin of the relatively poor quality of our hts, 
we have measured yet another observable, the average linear extension 129 


1 




r = — > r • n(r] 
N2 


(13) 


of a given triangulation T of discrete volume iV 2 , which is just the weighted aver¬ 
age of the geodesic distance r. As is described in more detail in the Appendix, the 
observable f has convergence issues, which appear to persist even on large lattices 
and after a large number of sweeps. What seems to happen to the geometrical 
conhgurations is that most of the time they are approximately “square-shaped”, 
with comparable linear extensions for either torus direction, but ever so often 
make an excursion to an overall shape that is elongated, where one torus di¬ 
rection becomes longer than the other one, with f increasing as a result. After 
a relative maximum of the two lengths has been reached, the system gradually 
reverts back to being square-shaped and stays there for a while before another 


excursion takes place (see Fig. 16 for illustration). 

Of course, square and elongated configurations (for identical volume N 2 ) not 
only have different average extensions f, but also different shape functions n(r) 
(see Appendix) and therefore in general different Hausdorff dimensions. A likely 
explanation for our inaccurate determination of dh is therefore the failure of the 
shape to stabilize during the course of the simulation, and the hnite-size effects 
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associated specifically with elongated shapes, in addition to the already mentioned 
large magnitude of the quantum fluctuations overall. This is supported by a 
numerical experiment we have performed in pure CDT quantum gravity at a 
volume iV 2 = 9.000. For f'roT = 80 time slices, the behaviour was “square-like”, in 
the sense that the Monte Carlo history of the average extension f did not have 
any peaks. However, when we shortened the time extension to Itot = 20, peaks 
similar to those depicted in Fig. [TO] appeared. 



Figure 17: A bubble has a ssf-triangle at either end and arbitrarily many stt- 
triangles in between. 


Lastly, in our search for ways to improve the convergence behaviour of LCDT 
quantum gravity, we investigated what happens when self-overlapping bubbles 
are not allowed to occur. We mentioned these structures briefly in Sec. above. 
A bubble is a contractible loop of spacelike links, which in its interior is decorated 
by timelike links only (Fig. [Tt] ). It always has two sst-triangles at its end points 
and consists of sff-triangles otherwise. When a bubble winds around a compact 
torus direction, it can touch itself again (“self-overlap”) along one or more space¬ 
like edges of its boundary (see Fig. 18 for a simple example). Self-overlapping 
bubbles are geometrically significant, because - depending on their interior ge¬ 
ometry - they can give rise to timelike cycles as dehned at the beginning of Sec. 
1^ Their appearance is not forbidden by local vertex causality. Relevant to our 



Figure 18: A self-overlapping bubble; the links with label a are to be identihed. 

present discussion is the fact that globally self-overlapping bubbles cause severe 
thermalization issues in 2 -|- 1 dimensions, and therefore were removed from the 
ensemble [131 Ell- This motivated us to remove self-overlapping bubbles from 
the LCDT ensemble in 1 -|- 1 dimensions too, and to check whether it makes a 
difference to the measurement of the Hausdorff dimension. 
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Figure 19: LCDT quantum gravity without self-overlapping bubbles: Fit for best 
overlap of the rescaled shapes {n{r ))as function of the rescaled distance 
x = r/^ for Hausdorff dimension = 


Detecting whether a self-overlapping bubble is created during the Monte Carlo 
simulation (and discarding the corresponding move) is nontrivial, since the prop¬ 
erty is nonlocal and requires a computationally expensive walk around the lattice 
(see [ 23 ] for details on implementation). For this reason we performed the nu¬ 
merical analysis on slightly smaller lattices of volume iV 2 < 60.000. Measurement 
of the average linear extension f in this setting, at N 2 = 50.000, still revealed a 
peak structure similar to that of standard LCDT with self-overlapping bubbles, 
providing evidence that this structure is not responsible for the observed instabil¬ 
ity. Proceeding like before to determine the Hausdorff dimension for this system, 
via hnite-size scaling to maximize the overlap in shape (see Fig. 19) yielded a 
Hausdorff dimension of dh = 3.10 ± 0.2. 

To summarize, we have pinpointed an instability of the system with regard to 
its global behaviour, due to occasional excursions to a globally elongated state, 
which can be observed by monitoring the geometry’s average linear extension f. 
This is the likely source of the suboptimal data quality for the measurement of 
the Hausdorff dimension dh- By considering a more general htting function for 
extracting dh and by using a modihed ensemble without self-overlapping bubbles 
we have found no hints of additional sources of error or a shift of the Hausdorff 
dimension toward the CDT value of 2. 




7 Conclusions and outlook 


The aim of our work was to measure observables in locally causal dynamical tri¬ 
angulations in two dimensions, most importantly, the spectral and the Hausdorff 
dimensions, and thereby understand the relation of LCDT to other models of 
two-dimensional quantum gravity based on dynamical triangulations. Our initial 
hypothesis was that LCDT lies in the same universality class as CDT, where 
both spectral and Hausdorff dimension are equal to 2. While our measurement of 
LCDT’s spectral dimension did yield a value compatible with 2, with only small 
error margin, this was not true for the Hausdorff dimension. Although the error 
bars were signihcantly larger - due to an instability in the system that persisted 
even at the largest volumes - our measurements found a Hausdorff dimension of 
dh = 2.7 ± 0.2 and dh = 3.1 ± 0.2 for two slightly different variants of LCDT. On 
the basis of our simulations, it appears that LCDT is not equivalent to either DT 
or CDT in the continuum limit. 

This would be an interesting result, because it implies the existence of a new 
universality class of two-dimensional quantum gravity in between Euclidean DT 
and Lorentzian CDT in two dimensions. The “in between” could be true quite 
literally, since within our measuring accuracy the Hausdorff dimension of LCDT 
is compatible with 3. A Hausdorff dimension dh = 3 has been observed previously, 
in simulations of CDT quantum gravity in 1-1-1 dimensions coupled to eight copies 
of Ising spins [32] and coupled to several massless scalar helds [33] , adding some 
plausibility to the possibility that a universality class with this property may 
actually exist. Further conhrmation of the appearance of this new phase would 
come from locating a phase transition between CDT (corresponding to a = 1/4, at 
least for hxed volume) and LCDT as a function of the parameter a in our model. 
Having already invested considerable computing resources into the isotropic case 
a = 1 in the present work, we leave this investigation to a future publication. - 
Needless to say, it would be extremely interesting to hnd an analytic solution of 
the LCDT model, to put our Endings on a more definite footing. 
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Appendix 


In this appendix, we give some more details about the instability we have ob¬ 
served in the LCDT system, which affects at least one observable, the average 
linear extension f of the universe defined in (13), and contributes to the rather 
poor overlaps we have found in our hnite-size scaling analysis to determine the 
Hausdorff dimension. Unlike other observables, which typically converge after 
about 300 sweeps, f does not, even for very large system size N 2 = 400.000 
and after several thousands of sweeps. To understand better what happens ge¬ 
ometrically, we have plotted the shape n{r) of a typical conhguration along the 
meta-stable “bottom” of the Monte Carlo history of f shown in Fig. 16, and of 
a conhguration at one of the peaks. As illustrated by Fig. the two are very 




Figure 20: Average shape n(r) (top) and shape at peaks in r (bottom) at A^ = 
9.000, exhibiting a long plateau. 


different. Outside the peaks in f, the shape of a conhguration starts out with 
an almost linear increase until it reaches a single maximum, and then quickly 
drops to zero. A conhguration from a peak in f also increases linearly until a hrst 
maximum, but then enters a long plateau before also going to zero. These two 
diherent shape functions are characteristic for a torus which is approximately 
“square-shaped” (i.e. of a similar extension in either of the torus directions) 
and one which is elongated. This fact is illustrated by comparing the measured 


shapes with those of regular, hat tor:^^ in the continuum (Fig. 21). Despite the 


^"^For simplicity, we are considering only tori which are obtained from gluing flat rectangles 
without any “twists”. 
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totally different set-up (single, classical torus without local curvature and without 
quantum fluctuations), there is a clear qualitative resemblance with the shapes 
extracted from the full quantum simulation. Note that we have not attempted a 
proper translation between discrete and continuum units of length and volume, 
which would be necessary for a quantitative comparison. 


n(r) 



n(r) 



Figure 21: Shape of a flat, classical torus in the continuum: square-shaped of 
length 50 in either direction (top) and elongated with extension 30 and 180 in 
the two directions (bottom). 
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